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SUMMARY 


The unsteady laminar boundary layer induced by the flow- Initiating 
shock wave passing over a flat plate mounted in a shock tube is theoretically 
and experimentally studied in terms of heat transfer rates to the plate 
for shock speeds ranging from 1.695 to 7.34 km/sec. The theory presented by 
Cook and Chapman for the shock- Induced unsteady boundary layer on a plate is 
reviewed with emphasis on unsteady heat transfer. A method of measuring 
time-dependent heat-transfer rates using thin-film heat-flux gages and an 
associated data reduction technique are outlined in detail. Particular con- 
sideration is given to heat-flux measurement in short-duration ionized shock- 
tube flows. Experimental unsteady plate heat transfer rates obtained in 
both air and nitrogen using thin-film heat- flux gages generally agree well 
with theoretical predictions However, some difficulties were encountered in 
resolving time dependent heat transfer rates for very early times after 
shock wave arrival at the gages. The experimental results indicate that the 
theory continues to predict the unsteady boundary layer behavior after the 
shock wave leaves the trailing edge of the plate even though the theory is 
strictly applicable only for the time interval in which the shock remains on 
the plate. In addition, the experimental results confirm that under certain 
conditions steady flow can readily be attained over models in short-duration 


shock-tube flows . 
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SYMBOLS 

Shock strength parameter, <jV/P ^fo^l) ' 1 > 13 U </ U D 

speed of sound 
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correction factor in Eq. (59) 

initial voltage drop across thin-film hear flux gage 

time dependent voltage change across thin-film heat flux gage 

enthalpy 

dimensionless enthalpy, (h - h 0 )/h Q 
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’ jr o o / o 

ratio. Fig. 4, defined by Eq. (55) 

ratio. Fig. 4, defined by Eq. (54). 

resistance of ballast resistor, Fig. 6 
change in thin film resistance 
initial resistance of thin film 
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time, t 1 — 0 v. 7 h eu shock wave arrives at thin film heat flux gage 

temperature change 9 - 9^; also absolute temperature when subscripted 
as T x 

velocity in the x direction 
shock speed 

velocity in the y direction 
width of thin film element. Fig* 6 
distance from plate leading edge 
distance measured from the shock, Fig. 2(c) 
gage position. Fig. 12 

distance perpendicular to plate surface 

distance perpendicular to wall, Fig. 2(c) 

dimensionless distance- time variable, x/u Q t 

thermal diffusivity of coating material 

thermal diffusivity of film 

dimensionless velocity, u/u Q 

dimensionless velocity expression, Eq. (14) 

1 n 

dimensionless time also a Reynolds number based on x - u Q t 

resistance- temperature parameter, AR /,A®R 0 

coating thickness, Fig. 6 

film thickness, Fig. 6 

general x in Eqs. (6), (7), and (8) 

absolute temperature 

C< F t/ > Fig. 7 
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film 
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INTRODUCTION 


Shock tubes have been developed that have the capability of driving 
f low- initiat ing shock waves at speeds in the neighborhood of 10 km/sec. 

The high- temperature high-speed flow behind the incident shock is of 
interest in a number of areas. Experimental studies in many of these 
areas are influenced by or are involved with the boundary layers induced by 
the shock wave on either the shock tube side wall or on models or test 
apparatus such as splitter plates placed in the flow. A difficulty 
encountered in shock tuba testing, however, is the short duration of the high 
tempei’ature flow. This is especially pronounced at high shock speeds. In 
view of the short testing times available, an important consideration is the 
transient development of the shock-induced boundary layers since, in cases 
where, steady-state boundary layers might be desired, the testing time may 
be insufficient for steady state conditions to be reached. 

An unsteady shock-induced boundary layer of considerable interest is 
that developed on a flat plate with a sharp leading edge mounted with zero 
angle of attack in a shock tube. The present study is concerned with this 
boundary layer. Considerations here are restricted to laminar boundary 
layer flow. The development of the laminar boundary layer on the plate has 
been the subject of several investigations, the most recent one being that 
of Cook and Chapman 1 in which a complete theoretical description of the 
unsteady boundary layer is presented for shock speeds ranging to 9 km/sec. 
However, experimental vertification of the theory is lacking for a 
large part of the range of shock speeds considered. The present investi- 
gation was initiated for the purpose of experimentally verifying the theory 
of Ref. 1. A quantity associated with the transient nature of the boundary 
layer that can be measured is the plate surface heat flux. As a result, 



?. 


this investigation focuses primarily on experimental verification of 
theoretically predicted heat transfer rates for the unsteady laminar flat plat 
boundary layer. The experiments reported here cover shock speeds ranging ' 
from 1,695 to 7.34 km/se.c and were performed using air and nitrogen as the 


test gas. 



theory of the unsteady shock 


INDUCED LAMINAR BOUNDARY LAYER ON A FLAT PLATE 


In view of the complexity of the theoretical analysis of the shock 

induced unsteady laminar boundary layer on a flat plate, the theory presented 

in previous papers on this subject will be summarized here in some detail. 

2 A 

Larn and Crocco originated the study of the shock induced boundary layer 

5 S 

on a flat plate. Feldennan 5 extended the work of Lam and Crocco by 
obtaining an extended range of theoretical solutions to the governing equations 
to predict unsteady heat transfer rates to the flat plate. He also experi- 
mentally verified the theory for a range of relatively lov/ speed flows. More 
recently the basic approach of Lam and Crocco x*as used by Cock and Chapman 1 


in the study of laminar boundary layers induced by shocks with speeds ranging 
from 1*12 to 9 km/sec, A complete description of the boundary layer develop- 
ment is presented for this range of shock speeds for equilibrium air at a 
pressure of one atmosphere behind the shock. Results describe the unsteady 
nature of the boundary layer in terms of friction coefficients, heat transfer 
quantities, and velocity, displacement, and momentum thicknesses as a 
function of a single position- time variable. The following sections present 
a review of the considerations that served as the basis for the analysis 
presented in Ref. 1. 
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Boundary Layer Description 

The boundary layer of interest is that shown schematically in Fig. 1(a). 

The shock wave is assumed to propagate into region 1 with constant velocity 
U s> At zero time the shock is at the leading edge of the stationary plate 
and by time t the shock wave has traveled a distance t in the positive x 
direction, and an unsteady boundary layer has developed on the plate between the 
leading edge and the shock. The free-stream velocity behind the shock relative 
to the plate is u q , and in time t the free-stream fluid particles initially 
over the plate leading edge have moved a distance x = u^t. The region between 
U s t and u o t contains only fluid particles that initially occupied positions 
downstream of the leading edge. Thus the boundary layer developed in this 
region (denoted as the >1 region) is identical to that induced by the moving 
shock on the shock tube side wall since the plate leading edge does not 
influence M region particles. Near the leading edge of the plate the 
boundary layer is that for steady flow over the plate. As will be seen, 
these two regions play an important role in the solution of the governing 
boundary layer equations. 

A key factor in the physical interpretation of the boundary layer develop- 
ment and the solution of the governing equations is selection of a coordinate 
system in which to view the problem. The variables of Lain and Crocco^’^ 
form a "natural" set of independent variables that readily permits interpre- 
tation of certain important features of the problem. First, a transformation 
is made to change from the variables of Fig. 1(a), (x,y,t), to (x,u,t) where 
u is 'the x component of velocity. Next, dimensionless variables are defined 

as b! = x/u 0 t, G= u/u , and \ ~ u t/T*. , where is the kinematic 

\ o o o o 

viscosity. An interpretation of the problem in terms of these variables is 
presented in Fig. 1(b). ft = 0 and ~ 1 correspond to the wall and 
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free- stream respectively. The region of interest on the o>; axis ranges 

from zero to A. U t/u t - U /u - A defines the shock position. The 

s o s o 

region 1 a, v /> 0, is the M region and » 1 delineates the dividing 

line in the free-stream that separates the fluid particles initially upstream 

r 

and downstream of the leading edge. The region 0*Zo( £ 1, 0£ ^ i 1, V>0, 

is denoted as the interaction region after LanTh Small values of correspond 
to steady state conditions. This can be seen by considering the problem 

from two points of view. For any fixed x position, increasing time corresponds 

j 

to decreasing values of and approach to the steady state condition for 
that x. corresponds to infinite time, Alternately, v?h en time is 

fixed and small values of x compared to u^t are considered (i.e., at small ^ )* 
steady state would be expected. Thus the dual role of the &< variable 
becomes evident. As noted schematically in Fig. 1(b), the region of interest 
in OC t ^ J V" space is 0^ ^ :£ 1 , 0 6 ot £ A, 0, Prior to considering 

the governing equations for the unsteady problem, solutions for the bounding 
M ana steady state regions will be discussed. 

Solution of the Steady State and M Region Boundary Layer Equations 
A unique similarity exists between the steady flat plate boundary layer 

and that for the >i region when, the M region is considered in shock fixed 

coordinates in that the governing equations and the boundary condi cions for 

these seemingly different cases are nearly identical. The governing equations for 

* 

for these boundary layers have been combined into a single set of equations in 
Ref. 7 which is’ summarized as follows. 

Consider the steady boundary layer on- a stationary flat plate as shown 
in Fig. 2(a). The wall velocity u is zero and the free- stream velocity is 
denoted by u f0 , A typical velocity profile is shown. Fig. 2(b) shows a 
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shock wave moving with velocity U into region 1 along and perpendicular 

v> 

to a flat wall of infinite extent. The velocity of the fluid behind the 
shock is U2- The shock wave may be fixed by subtracting U g from every part 
of the flow field, including the wall. The resulting flow is illustrated 
in Figure 2(c). Note that the wall velocity u w is greater than the free- 
stream velocity uoo • Also, the transformation to shock-fixed coordinates 
removes the unsteadiness which is present in the wall-fixed reference frame. 
The similarity between the M region and steady-state boundary' layers is 
now evident. The only difference is that the wall velocity is finite for 
the M region case and zero for the steady-state case. 

The analysis may continue now by treating the steady-state boundary 
layer. Fig. 2(a), and the M region boundary layer, Fig. 2(c), as the same 
problem except for the wall velocities. The continuity, momentum, and energy 
equations for steady, thin two-dimensional, compressible, laminar boundary- 
layer in thermochemical equilibrium with zero pressure gradient are: 


Continuity: 

h X 6 $ 

Momentum 


0 U- Aifr + pv A 'A- _ r ,, chx \ 
r ^ ^ 

Energy 

r ^ r h V 


a) 


( 2 ) 


(3) 


In these equations x is either x in Fig. 2(a) or x 1 in Fig. 2(c). In Eq. (3) 

8 

k is the total thermal conductivity . 
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The boundary conditions, referring tc Figs, 2(a) and 2(c), are: 


0 for steady- state case 


u(x,0) 

u VJ for M region case 

(4 a) 

v(x,0) 

- 0 

(4b) 

h(x,0) 

N-? 

(4c) 

h(x,os ) 

« 

(4d) 

u(x,o> ) 

= 

(4e) 


In Eq. (4c) h w is assumed constant. Using boundary conditions (4a) and 
(4b) ? evaluation of Eq. (2) at the wall (y = 0) yields 
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O 

where cr is the effective Prandtl number^ which implicitly incorporates 
the influence of dissociation and ionization. The transformed boundary 
conditions are: 


at u = u w (y “ °) 

h = h 

w 

V* = 0 

u 

at u = (y =<x>) 

h = h„ 

CO 

o 


The following dimensionless variables are now introduced: 


r * 


r 




pco U“<*> 


? S 


u w - <£ 

U«, 


H = 


h-h 


CO 


l CO 


The momentum equation then becomes 


(9a) 

(9b) 

(9c) 

(9d) 

( 10 ) 

(ID 

( 12 ) 


! y -r — ) - o 

(I 


( 13 ) 
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where 



Note that for the steady state case 



since u * 0. 
w 


The energy equation becomes 


( 14 ) 





The boundary conditions become 

at ^ 0 and all ^ 

^ = 0 
H = H w 

at ^ - 1 and all ^ 

r - 0 

H - 0 


(16a) 

(16b) 


(16c) 

(16d) 


Note that the difference between the M region and the steady state 
boundary layers is incorporated in the equations in terms of i^. The 
boundary conditions are now the same for both cases. 

The solution of the momentum equation is found by separation of variables. 
The solution for the shear stress is: 
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where 0 (^) is the solution of: 

M., = O 


W 


Z 


and c » 


Substituting Eq. (17) into Eq. (15), the energy equation become: 


( 18 ) 

(19) 



j2 (u.co “ Llw) 




( 20 ) 


i 

It is noted that the form of Eq. (20) does not change with increasing j 

: w» 

and that the boundary conditions on H, Eqs. (16b) and (16d) , are the same 
for all ^ . Thus the solution does not depend on ^ j and the last term 
may be dropped from Eq. (20). The energy equation can then be written as: 



The boundary conditions are now 
- 0, H - H w at $ « 0, 

0 •** 0, H = 0 at 0 = 1. 

Thus Eqs. (18) and (21) may be written in integral form as: 


( 21 ) 




11 


H(p a - Hg('o> I 

p crfo) 1 


+ 


(5 


(U.a>-U w f f (J-CP) rP 


• C3 




5<W 


§(N)dHd'P 


+ 


cr(P> 




(p (p) \ ^ [$(**)] CT^ cfA/d P 

^ ‘ Jo 


(23) 


where 





l-0-(?> 


and P and N are dummy variables for ^ . Coupling exists between Eqs. (22) 
and (23) since C, Eq, (19), appears in Eq. (22) and is a function of H. 

Of particular interest here is the wall heat flux q w » -k( d T/ ^ y) , 

which, when written in terras of the present variables, becomes: 




‘T'*/ A « 

o; v (u m -u w ) v ‘fV 


(H s l 


(24) 


where (H ^ is the enthalpy gradient at the wall ( ^? =* 0) . Substitutes 
for K from Eq. (17) a dimensionless heat transfer quantity is obtained. 



( 25 ) 
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/ 


Thus, specification of in Eq. (22) (through Eq. (14) ) and in Eq. (23) 
permits evaluation of the dimensionless shear stress, F.q. (17), and the 
dimensionless heat transfer, Eq, (25), for either the steady state boundary 
layer or the M region boundary layer viewed in shock-fixed coordinates. 

At this point interpretation of the quantities obtained above can be 
made in term of the coordinates in Fig. 2. Consider first the steady 
boundary layer in Fig. 2(a). The profiles 0(^) and H(^) obtained from 
Eqs. (22) and (23) with u w » 0, denoted here as 0(^) ss and H(^) ss , apply 
to this boundary layer. In addition, the shear stress and heat transfer 
from Eqs. (17) and (25) are 



(26) 



V3T <r w 
h 0 \j p (i yiolio 



(27) 


where ^ in. Eqs. (17) and (25) is replaced by x, the plate distance 
coordinate in both Figs. 1(a) and 2(a), and the subscript n o H is used to 
denote the free- stream quantities. 

Next consider the M region boundary viewed in shock fixed coordinates, 
Fig. 2(c). The 0(^) and H(^) profiles obtained from Eqs. (22) and (23) 
with u^ =* which are denoted as and a Pply to the 

boundary layer in this figure. The expressions for shear stress and heat 
transfer from Eqs. (17) and (25) are: 


r(v:,'(3) - 


( _ Ik \lithi l — T ' d>($) 

V UaJSJ> a Ut,X ' Ms# 


( 28 ) 
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Q. 


Ms-? 


^ w /~ X 0 \*j 

CO \l J~bo jJL 



( 29 ) 


where x* is the distance coordinate measured in the downstream direction 

relative to the shock as shown in Fig. 2(c). It is also noted from Fig. 2(c) 

that tWu^ - the density ratio across the shock. 

When a point on the moving wall is at the origin of the x' , y' coordinate 

system of Fig. 2(c) at time t' = 0, the position of the point at time t' is 

x' =■ U t'. Note that this is also the distance the moving shock travels 
s 

past a point on the shock tube sidewall in time t'. Substitution of x' - 
into Eq. (29) yields 


o\/7’= 


' w 


a 




J Msf 


(30) 


given 


Thus, for a ^sidewall boundary layer, q \fF = const. 

Eqs. (26) and (27) and Eqs, (28) and (29) provide the shear stress and 
heat transfer solutions in the coordinate systems noted for these equations. 
It is next of interest to transfer each of these solutions to the coordinate 
system of Fig. 1(b), i.e,, to the coordinates o(~ x/u^t, ^ «*u/u q , and 
Y “ u^t/y^. For the steady state region this is readily accomplished and 
Eqs, (26) and (27) become: 


=■]— nr 


(31) 


Qj 5 ,0, \ A 1 r 


o w y i 




c Kf 


f» u< 


! A 

- iV 


w 


w 


(32) 
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For the M region of Fig. 1(b) it is noted that x + x' = U t, where x is 

s 

the distance coordinate in. Fig. 1(a) and x' is the distance coordinate in 

Fig. 2(c) in which u is the free stream velocity. Hence u + u = U . 

oo J o co s 

Substitution for x 1 and u^ in Eqs. (28) and (29) yields for the M region: 



(33) 


(34) 


Further interpretation of Eqs. (31) to (34) will be made after the 
interaction region solutions are discussed. 
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Solution of the Boundary Layer Equations for the Interaction K^siotj 

The equations for the lroninar boundary layer in the interaction region 

i 

for real gas flows in th ervr.o c h era ic a 1 equilibrium are: 


continuity 


momentum 


A 0 

bt 

c X.p!X) . 


d ; >( * 

6 M 

pi a 

* dt 

i s*> if o lx 

+ r u gy 





~ O 


<^LL 


/> V\ 


energy 


p £b. -i- p u, ii + pit 4h 

J&t i dvc *)<f 


agt ^ 


(35) 


(36) 


(37) 


where k in Eq, (37) is the total thermal conductivity and the coordinate 
system is that of Fig. 1(a). From previous discussion it is noted that 
the following boundary conditions apply. 

u(x,0,t) » 0 
v(x,0,t) » 0 
h(x,0,t) ~ h w = const. 
u(x,oo ,t) = u Q 
h(x,co ,t) = h D 

Evaluation of the momentum equation at the wall yields: 


(38a) 
(38 b) 
(36c) 
(38d) 
(38e) 


4- (u. | - ~- L Is o 

h <4 J \ 

o * u*o * '^*0 


oil 


- kfu 


( 39 ) 


tjS'G v *g*0 

which replaces boundary condition (38b). Transforming the equations from 
the independent variables (x,v,t) to ( % « x, u, £ * t) and incorporating 
the continuity equations into the momentum and energy equations, Eqs. (35) 
to (37) become: 
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momentum 


a (£*) + rL + f£&\ = o 

i >■ ; i f y • -» 

' i OU \ » / 


(40) 


energy 


/>P»? + a, yj +^*w 


(41) 


O'" in Eq. (41) is the effective Prana cl. number . Equations (40) and (41) written 
in terms of the independent variables of Fig. 1(b), 6( = x/u c t , - u/u q 

m 9 2 

and ^ =* u t / , with H - (h - h ) /h and 't* ** T / & }J as independent 

o * o 5 . o / o # J ° ° 

variables are: 


r\ - c(^jz + c v - rc r + 


r c « 


(42) 


c [*v - Ht (^)J + r f n f = !jj r% r(J 


(«> 


, A solution for the shear stress in terms of a new 


where C - 

shear stress variable 0 is defined as 


$y) = ^ 


(44) 


\fnen this expression is substituted into Eqs. (42) and (43), these equations 
become : 

+ Cj wfz = «r(c$r-$c,) (45) 


r u ' 
r-:3 


d AtiJ±~i) - 

i A> 3 7 


e,v.((3-°0 H + u|_ 




Cofy w 
“f" H v 


(46) 
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The region of interest: for the above equations is the interaction region 
0 £ « £ 1, 0 £ 1 , Y> 0, in Fig. 1(b). Hence, boundary conditions must 

-be supplied for this region. The boundary conditions on 0 at Ow ~ 0 are 
furnished by the previously discussed steady state solution of the momentum 
equation. Lam~ has shown that the solution at o( = 1 is that for the M 
region at c< ** 1. Combination of Eqs. (33) and (44) yields 

which expresses 0 for the M region of Fig. 1(b). The boundary conditions 
for 0 and ^ = 1 are obtained by transferring those listed in Eqs. (38) 

to the (c*',R , Y) coordinate system. The boundary conditions for Eqs. (45) 


and (46) then become: 



0 (0, |3,Y) = 0(f ) ss 

H(0, V) - H(e) gs 

(48a) 

0 d.p.Y) -flKp) M8f 

H(l, ) = * l ( ^ ) ]q s f 

(48b) 

/ 

0 («,o,y ) - K 

H(trf,0, V) = H w 

(48c) 

0 (o< ,1, Y ) = 0 

H(o< , 1 , Y ) « 0 

(48d) 

It is noted that these boundary 

conditions are independent or * . 

Further, 


it is noted that Eqs. (45) and (46) do not change in form as Y changes, i.e., 
when Y* ~ B Y (where B is an arbitrary constant) is substituted into Eqs. 
(45) and (46), the expressions on the righthand* side of these equations 
undergo no change in form. Hence 0 and H are independent of Y and therefore 
the » derivatives are zero. H and 0 remain as functions of and ^ only. 
Integration of Eqs. (45) and (46) along lines of constant e< along with the 
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application of the boundai*y conditions in Eqs. (48) yields: 


r rJ 




J 




*6M> 1 


A - 2 *S( A - c<) 


/* 

' 


C« (W aV 
C^A) 


>'A .*? 

ex A w) 


j 




(h?) 


H (K?) 


H 6<,o) ± / JgJgLil 

> ’ ' crfrs) J *l- 


4f 


*' "kf) 




) | C<& M {A ~ ^ ) Ho. { J v ^ ^ 

\i + <r(«,x\. ; 

A C^h) 


,> oTe^i 

♦ <W,$) 


'o 


JA o I or*^/.;) 




'e 


A l -&(<*>% I . 

<t> <r*,*) 


4> 


| - G*(U ; A> 




v Q 


+ 




q~iX^ 

$ 'V^) 


cp A <*/) a A 4 f 


cnf*/) 


(50) 


"^ef, 1 lists four additional ’’derived" boundary’ conditions which aid in the 
numerical solutions of Eqs. (49) and (50). They are: (j) -- O ^ 
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As was the*, case for the energy and momentum equations for the steady state 
and M regions, Eqs. (49) and (50) are coupled since. C in Eq. (49) is a 
function of H determined by Eq. (50) and 0 in Eq. (49) also appears in Eq. 

( 50 ) . - 

Of concern here is the expression for the wall heat flux q^ = -k(^T/ <i y) w 
for the interaction region. In terms of the present variables, 


U M r 


- u< 

Oh/ 




(51a) 
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(51b) 


where 0(e<,O) and (t<,0) are from the solutions of Eqs. (49) and (50) 

respectively and substitution for ? * from Eq. (44) has been made. Equa- 

' 1 

tion (51b) can be written in dimensionless form as 


Q (o<, o) 



A u » 




(52) 


Thus the heat transfer variable Q has been specified for the complete range 
of interest: referring to Fig. 1(b) , Eq* (32) applies to the steady state 

portion of the interaction region, Eq. (52) is applicable to the entire 
interaction region, and Eq. (34) applies to the M region. It is noted that 
Q depends only on the time- posit ion variable o4. 
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Theoretical Heat Transfer Solutions 

In order to obtain the desired solutions for a given flow condition, i t 

is necessary to numerically integrate Eqs. (22) and (23) for the steady 

state and M regions and Eqs. (49) and (50) for the interaction region. The 

method of solution of these equations is outlined in Ref. 1 and 7 and need 

not be discussed. Of the several descriptive boundary layer characteristics, 

only the solutions for the heat transfer quantities will be presented^ since 

heat transfer to the plate is the primary concern of the present study. 

When real gases in equilibrium at high temperatures are considered, 

the dimensionless solutions to the boundary layer equations depends on pressure, 

since pressure influences the variation of the product ratio C * 

the effective Prandtl number <T , and the density ratio across the 

incident shock. The solutions presented in Ref. 1 were obtained for air at 

a free- stream pressure p Q of 1 atm and T w - 300°K, The relations for viscosity 

and density variation with enthalpy used in obtaining those solutions are 

9 10 

those presented in equation form by Mirels and Viegas and Howe . Mire! s' 

equations were used for values of h/h up to 4.53. The equations of Viegas 

w 

ana Howe, which are based on Ref. 8, were applied at. higher values of h/h w . - 
The effective Prandtl number variation with enthalpy was taken from Ref. S. 
Free-stream properties behind the incident shock were obtained using Refs. 11, 

12, and 13. The solutions for heat transfer % at P Q ~ 1 atm from Ref. 1 serve 
as a basis for the discussion in the following paragraphs. 

As noted in Ref. 1, interpretation of the solutions in terms of the 
present variables for the boundary layer 'Considered here can be made by 
taking two different points of view. Interpretation for fixed x positions 
is desired for the present study because measurement of the heat flux to the 
plate is made at fixed x values. Also of some interest is the nature and 
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configuration of the boundary layer at any point in time. However, for the 
reason noted, only -the fixed x interpretation will be made here. It is 
noted that the quantity (c4V)^ appears in each of the equations, Eqs. (32), 

(34), and (52), that collectively predict the heat transfer quantity Q vs oC. . 

* 

, 1 

Substituting for oi and (ocV)'* - (u 0 x / 2 =* (Re)^, Thus for fixed x, 


= gw (ZzY^/<y'oU*K) - 


?.v/Vx orv 


\Jpo JX < 




( 53 ) 


where Re is the Reynolds number and x is the distance measured from the 
plate leading edge. A typical solution for Q x vs e< is presented in Fig. 3. 

When x is fixed, Oi - x/u Q t becomes a time variable and ctf decreases with 
increasing time. The time x/U s corresponds to the shock arrival at the 
given x location and the corresponding value of is A. Times less than 
x/U s (values of oi > A) are of no interest. With increasing time the boundary 
layer develops and approaches the steady state condition for the given x 
location. This behavior is observed in terms of Q x in the figure. Eq. (34) 
for the M region predicts the expected infinite value for Q x at ci = A. With 
decreasing ex. (increasing time) Q x decays to (0 W 0Hft) w )>i S f at ot = 1 according 
to Eq. (34). Q x continues to decay along the curve predicted by Eq. (52) and 
approaches the steady state solution, Eq. (32), asymtotically with decreasing oi . 

Several important observations can be made regarding the solutions for 
Q x typified by Fig. 3. First, a given solution Q x vs ot is self similar -in 
the sense that it applies at any fixed x provided the plate is at least of 
length U s t. Next, as might be expected, the steady state Q } . is essentially 
reached at values of c*' > 0, i.e., before infinite time. Further, the 
interaction region solution is asyvnototic to both the steady state and M 
region solutions as is illustrated by the extended curves for these regions. 
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Use has been made of this fact to correlate, the 

Fig. 4 presents such a correlation for air. 

defined as „ /c „. 

R * 

ss Q >: , Eq. (32) 


v> = Qx> E * < 52 ) 

™ Qx, Sq. (34) 


beat transfer results, (Kef. l). 
This figure presents two ratios 

(54) 

4 

(55) 


in terms of and the chock speed U s . This figure permits Q x vs to be 
predicted provided the steady state and M region values for Q x are available. 
Quantities which permit the determination of Q x for the steady state and M region 
for air at a pressure of 1 atm are presented in Fig, 5. To construct a 
solution Q x vs o< for a given U g , the known steady state and M solutions are 
extended in the manner shown in Fig, 3. Values of the ratios R ss and Rj.j 
and the corresponding values are then determined from Fig. 4 at the given 
U s , Eqs. (54) and (55) are next used to predict values Q x (at the several 
os. values corresponding to the ratio values in Fig, 4) for the interaction 
region. The interaction region curve can then be drawn through these points 
asymtotic to the steady state and M region extended curves. This method 
eliminates the somewhat tedious and expensive procedure of solving Eqs. (49) 
and (50) in order to compute Qx for the interaction region by Eq. (52), 


The inflvience of pressure on values for the dimensionless boundary layer 
quantities is most pronounced for shock tube flows generated by incident 
shocks at high speeds. For air with T w « Tj. - 300°K, the influence of pressure 
on the dimensionless boundary layer quantities begins at about U g = 2 km/sec, 
and increases with increasing shock speed. Therefore, for cases for which 
the incident shock speed is low, the solutions Q x vs c4 determined using Figs. 

4 and 5 would be satisfactory regardless of the value of p . 


For high shock 
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speeds where differs significantly from 1 atm., the quantities in Fig. 5, 

especially those for the M region, are not satisfactory, and it is therefore 

necessary to obtain solutions to Eqs. (22) and (23) at the given pressure 

that lead to computation of Q for the steady state and M regions by Eqs. (32) 

and (34) . For exact solutions in the interaction region at high shock speeds 

it is also necessary to solve Eqs. (49) and (50) at the given pressure in 

order to determine Q in the interaction region using Eq. (52). However, 

X 

since it is not easy to obtain these solutions, a good approximation of the 
solution for in the interaction region can be obtained using the ratios 
in Fig. 4 provided the quantities for Eqs, (32) and (34) are obtained at the 
given pressure. 

Fig. 4 is strictly applicable only to air. However, it is doubtful that 
the ratios R flg and would be significantly different for nitrogen. 
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HEAT FLUX MEASUREMENT IN SHOCK TUBE FLOWS 

In order to measure transient heat transfer rates in shock tube flows 

it is necessary to use an instrument with very fast response; one with a 

response time of a few microseconds or less. This consideration is particularly 

important when the flow duration is of the order of tens of microseconds. In 

addition to fast response it is important that the instrument be reliable, easy 

14 

to construct and use, and durable. The thin-film heat flux gage is the only 
practical type heat flux gage that satisfies most of these requirements and 
was therefore chosen to measure heat transfer rates in this study. The central 
idea of the gage system is to measure the temperature on or near the surface of 
a semi-infinite solid initially at uniform temperature exposed to the shock 
tube flow and to deduce the heat flux to the solid using heat conduction 
theory and the measured temperature. The purpose of this section is to des- 
cribe some of the important theoretical and practical considerations involved 
in heat flux measurement using thin-film heat flux gages. 

The Thin-film Heat Flux Gage 

Fig, 6 shows a schematic of the thin-film heat flux gage. A platinum 
film of thickness 6p is deposited on the surface of a substrate that is of 
sufficient thickness to behave as a semi-infinite medium during the shock 
tube testing time. Code 7740 Pyrex is a commonly used substrate since its 
thermal properties are well known. A coating'of thickness <5 C (usually silicon 
dioxide^) is deposited over the platinum film to eliminate electrical 
shorting of the gage if the gage is to be used in ionized flows. The platinum 
film serves as a resistance thermometer and is placed in a "ballast" circuit 
in which the electrical current is maintained essentially constant by making 
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very large compared to the film resistance* An oscilloscope and camera 
is used to record the film voltage change vs time during testing* The gage 
voltage response may be related to temperature change as follows. Considering 
the current 1 Q in Fig. 6 to he constant, the voltage change across the film 
is given by 

A E P = I q AR f (56 ' 

where ARp is the change in film resistance which is related to temperature 

change X to a high degree of accuracy in the range of temperature change of 

i 

concern here by the linear relation 
£» R p = (AR/A0)T 

where AR/AG is the slope of the resistance vs temperature plot for the 
film* Substituting into Eq* (56H 

AE F - I Q (AR/A9) T - I 0 R 0 (ARM0R o )T - E q P T 
or 

T * AE f /E 0 r (37) 

E q is the initial voltage drop across the film which has an initial 
resistance R q . Thus when I Q is constant* the temperature change T indicated 
by the film is related to E 0 , the film voltage change AEp (a function of 
time) recorded by the oscilloscope , and the resistance- temperature parameter P 
The latter quantity is usually determined by static measurement of resistance 
vs temperature for the film. 


*A more detailed analysis of the thin- film gage circuit which does not assume 
constant current yields the following expression for the temperature change 

T = [ 1 + <R 0 +AR p )/R 3 } <-AE f /E 0 r ) 

Thus if Rj}>100 R 0? Eq. (57) yields temperature change values that are 
correct within one percent. 
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Interpretation of Gage Response 

In general it is a complex matter to interpret thin-f ilm. gage response 
to determine the surface heat flux, especially for times of the order of 
microseconds. This is due in part to the finite thickness of both the 
coating and the platinum film. It is noted that the film would, due to the 
presence of the coating, respond to a temperature change that is not chat on 
the surface. A temperature gradient might exist through the film adding a 
further complication. No solution is readily available for this complex heat 
conduction problem. However certain special cases have been solved that permit 
insight into the influence of the factors mentioned above. 

It should be noted that the characteristic heat diffusion depth (see 
Carslaw and Jaeger, Ref. 16) is much less than the gage width W shown in 
section A- A of Fig. 6 for the time intervals of concern here. Therefore heat 
conduction in the gage materials can reasonably be treated in one dimension. 

Unco ated platinum film . - One special case of interest is that in which 
the coating thickness $ Q is zero, i.e., the case of an uncoated platinum film 
on pyrex. As a result of the fact that platinum thermal properties differ 
significantly from those of pyrex, the platinum film on the pyrex substrate 
will at very early times after heat flux application have a temperature different 
than the surface of the adjacent substrate not covered by the film. Hence at 
early times the film does not indicate the true substrate surface temperature 
denoted here as Kurzrock 17 and Camac and Feinberg 18 have analyzed this 

problem for certain cases of surface heat flux variation with time and 
determined the temperature change T at the interface of the platinum and the 
substrate. Fig. 7, which is based on the 'work of these investigators, presents 
curves Tj/T^ vs the Fourier modulus 6 for two cases of surface heat flux. 

The case of the surface heat flux proportional to t is of particular interest 
here because this is the heat flux variation with time that would be expected 
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on the shock tube side wall as the shock wave passes any given location. This 
heat flux case has not been investigated as thoroughly as has the constant 
heat flux case. However, it is noted that the approximate curves for this case, 
which are based on computations for a small range of 0p in Ref. 18, are not 
significantly different from those for the constant heat flux case. It is 
noted in Fig. 7 that in addition to a dependence on 9 = o/ t/ and a mild 
dependence on the nature of the surface heat flux, T^/T^ also depends on the 

F . (The subscripts s and F indicate substrate 
and film respectively and 0(p is the film thermal dif fusivity. In this section 
t is time measured from the instant of application of the heat flux and is not 
the time in the distance- time variable 0 ( .) For platinum on pyrex, = 0.1, 
as noted in the figure. The curves show that when 9 is very large (i.e. } for 
large times or small £ ) T /T is essentially unity regardless of the nature 
of the surface heat flux or the value of 4^ . Under this condition 

substitution of gage response (voltage vs time) into • - ' • 

/ 

Eq. (57) yields the true substrate surface temperature change. When substrate 
properties are considered to be constant and one dimensional heat conduction 
is assumed, the heat flux can be found by the expression^’ ^ 


quantity r = \fi kp7) s /(k j>c) 



j . 


(58) 


where T(0,t) is the substrate surface temperature change that would exist 

in the absence of the film (or coating) and is an arbitrary function of time* 

For the case where the modulus 0 is small, i.e., small times or a 

F 

relatively thick film, it can be seen in Fig* 7 that T^/T^ is less than unity. 
Fig. 8 presents for = 0.1 and the two cases of surface heat flux the ratio 
T T /T m vs time computed for a range of 0 for film thicknesses of 250 A° ana 

i- vV £ 
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1000 A°. Due to the high thermal dlffusivity of platinum, only a small- 
temperature gradient would bo expected through the film regardless of the value 
of 6 . Therefore the temperature sensed by the film would be essentially T^. 

This is then the temperature vs time that would be computed from the gage 

< 

response (voltage vs time) using Eq. (57). Since T is not the true substrate 
surface temperature, integration of Eq. (58j using this computed temperature 
would yield an incorrect heat flux variation with time. However, Fig. 8 can 
be used to estimate the true substrate, temperature as follows: 


T w - T(0,t) 


T(Eq. 57) 
cf 


where cf is a correction factor approximated by 


cf 




(59) 


( 60 ) 


and T /T is from a curve like those in Fig. 8 with ~ 0.1 plotted for the 
I co 

given film thickness and the approximate expected surface heat flux. Selection 


of a curve corresponding to the expected surface heat flux is necessary because 
the true heat flux is obviously not known a_ priori . Use of such a curve amounts 
to assuming that the nature of the surface heat flux is of secondary importance 
in determining T(0, t) . Support for the latter assumption is found in Fig. 8. 


It can be seen for the constant heat flux case that changing the film thickness 
by a factor of four has a pronounced effect on T^/T^ . However the influence 
of the nature of the surface heat flux is smaller, as can be seen from 
the curves for ^=0.1 for the two heat flux cases at a film thickness of 
1000 A°. Thus it is seen that precise knowledge of the film thickness is 
more important in determining cf in Eq. (59) by Eq. (60) than is kiiov?ledge of 
the nature of the surface heat flux. 

Coatin g over a very thin platinum film . Another case of interest is 
that in which the film thickness approaches zero but the coating thickness 
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is not zero. This would correspond to a coating over a very thin film. 

In this case the film senses the temperature a distance 6^ below the 
surface, not the true substrate surface temperature. Curves in Fig. 7 
permit analysis of this case also. Assuming that the coating has the sano 
thermal properties as pyrex, the modulus V — ~p c )j / c ) p ~ 1. The 
interface temperature change Tj is then given as a function of the modulus 
§ F (Qp = °<Qt/ Iq 2 in this case) by the '? - 1 curves in Fig. 7. Fig. 8 also 
presents Tj/T^ vs time for coating thicknesses of 250 and 1000 A 0 
for the two cases of surface heat flux. Since the film senses the interface 
temperature, T(0,t) for the case of a coating over a very thin film can he 
computed from Eq. (59) in which c f would be Tj/T^ obtained from curves like those 
in Fig. 8 with = 1, plotted for the given coating thickness and the approxi- 
mate surface heat flux variation with time. The curves in Fig. 8 for " 1 
again indicate that knowledge of the coating thickness is more important than 
precise knowledge of nature of the surface heat flux. After obtaining T(0,fc) 
in this manner the heat flux can then be evaluated by Eq. (58) . 

Determination of Heat Flux from Surface Temperature Variation 
Given T(0,t), the semi-infinite solid surface temperature variation with 
time, it is next of concern to determine the surface heat flux q(0,t) by inte- 
gration of Eq. (58), the integration of which is complicated by the fact that 
T(0,t) cannot usually be described by a simple mathematical expression. As a 
result it is necessary to resort to a numerical integration method. A method 
for numerically integrating Eq. (58) has been developed (Refs. 19-21) that yields 
accurate results with a minimum number of divisions of the time axis. This 
method approximates T(0,t) by a piecewise linear function which is introduced 
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into Eq. (58), When integration is performed the following expression is 
obtained: 


(oU , y°K- T(o,t,) - Tto.y,) 1 


1=1 (t la - t T r + (t m - 1 


m 


1-1 


—r; 

>* 


( 61 ) 


where tj is the time at the end of the 1th of m time intervals which are net 
necessarily equal. The thermal parameter \f(,kpz) i here applies to the 
substrate material. Due to diffusion of the film material into the substrate 
during the film deposition process the thermal properties of the substrate 
may be significantly changed. (See Refs. 22 to 24 for example.) Bogdan"^ 
investigated the variation o f for platinum films on code 7740 pyrex 

O i 

and determined the quantity to be 0,0737 Btu/ft °F sec 2 with a standard 
deviation from the mean of 0.0012. It is concluded that the mean value is 
sufficiently accurate for most purposes. Thus when T(0,t) is determined, 
discrete values T(0,tj) at time tj can be determined and introduced into Eq. 

(61) to numerically determine the heat flux. 

Influence of various factors on heat flux evaluation . - In order 

to illustrate the influence of various factors on the determination 

of the surface heat flux it is convenient to assume a gage response curve. 

Fig, 9(a) presents such a curve which represents a function that has an 

♦ 

asymtotic rise to a unit change in temperature. 

Time interval size: The first factor to be considered is the influence of 

the time interval size on the accuracy of the heat flux computed using Eq. (61) 
For this consideration it is assumed that the coating and film thicknesses 
approach zero. Thus the gage temperature response curve is the true surface 
temperature change T(0,t). The symbols in Fig. 9(a) indicate heat flux 
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values computed from the assumed temperature response using Eq. (61) and 
the time increments sizes noted in the figure for a period of about 21 usee. 

The heat flux is presented in terms of the quantity q^/t / |/(kpcJJ since it is 
of interest to compare the heat flux for the assumed gage response curve 
with that for a true step in surface temperature. Exact integration of 
Eq. (58) for a true unit step in surface temperature yeilds q|£ /y(Kpc) s - 1/ir 1 / 2 , 
which is noted on the figure. It is observed that the computed heat transfer 
values vs time are essentially the same for computations made using the two 
smaller time increments, indicating that a sufficient number of points have 
been used to accurately compute the heat flux for the assumed curve. At 
early times the computation of heat flux values using the 2.1 psec time in- 
crement differ considerably from those computed using smaller time increments. 

By about 20 usee results for all three time increment sizes agree well, and 
approach the theoretical value for a unit step in temperature as time increases, 
i.e., as the influence of the asymtotic rise of the assumed gage response curve 
diminishes. The results in Fig. 9(a) indicate that with sufficiently small 
time intervals accurate heat flux values can be computed using Eq. (61). 

An additional indication of the number of points required to compute heat 
flux within a prescribed percentage using Eq. (61) can be obtained by examining 
Fig. 10 which is taken from Ref. 19. This figure presents for the two boundary 
conditions noted the percent difference between the heat flux computed using 
Eq. (61) and the corresponding exact heat flux (determined from exact integration 
of Eq. (58)) as a function of the number of equal time intervals. It is noted that 
the boundary condition has an influence on the accuracy of Eq. (61). It is noted 
in Ref. 19 that the error in computing heat flux using Eq. (61) can be reduced 
by using unequal time intervals in such a manner that the boundary condition 
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T(0,t) is accurately described by a piecewise linear function. Thus in 
practice it is best to imagine a given T(0,t) curve as being approximated by 
a piecewise function and to choose the discrete values of T(0,t^) and t^. in 
such a manner that the curve is described as accurately as possible by 
straight-line segments. The errors for equal time intervals in Fig. 10 can 
be viewed as the approximate upper bounds for the respective cases noted 
therein that would be encountered using Eq. (61) . 

Finite film thickness; Another influence of interest is that of a finite 

j 

thickness platinum film. For this case, only the condition & = 0, as 

v 

j 

previously analyzed, will be considered* A film thickness of 500 A° is 

\ 

assumed. This value represents the maximum thickness expected for the films 

used in the experimental phase of this investigation and was estimated from 

film thickness measurements of platinum films on pyrex made by KcCaa, ref. 22. 

Due to the finite film thickness, the assumed gage response in Fig. 9(a) must 

be introduced into Eq. (59) in order to obtain T(0,t) and, in turn, discrete- 

values T(0,tj) and t^ for computation of heat flux by Eq. (61). The factor 

cf = Tj/T^ in Eq. (59) was obtained from a curve Tj/T^ vs t (like those in 

Fig. 8) plotted for p *= 0.1, - 500 A°, and qtt)^ = const. The ratio 

of the heat flux obtained for a film thickness of 500 A° (denoted as q(0,t) c ) 

to that obtained with no correction for film thickness (denoted as q(0,t) ) 

S r 

was computed for the assumed gage response curve in Fig. 9(a). This ratio is 
shown in Fig. 9(b). It is noted that the difference between the two heat flux 
values is not large except at early times and that the difference is less 
than 2% after about S^isec. v. 

Finite coating thickness: Another influence on heat flux evaluation from 
gage response is that of a coating over a platinum film. Again only the case 
previously considered, that of a coating over a very thin film, will be 
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analyzed. In this case the gage response curve. Fig. 9(a), is that for a 
very thin film placed between the coating of thickness S c and the substrate. 

A coating thickness of 1000 A 0 is assumed, and as previously noted, the surface 
temperature T(0,t) can be determined from Eq. (59) where cf for this case is 
approximated by the radio Tj/T^ for the 1000 A 0 coating, ip = 1, q(t) 1 ^ 2 = const in 
Fig. 8. Heat flux values for the corrected and uncorrected gage response 
curve, Fig. 9(a), were computed. The curve for the ratio of these two 
quantities is presented in Fig. 9(b). It is noted that the heat flux ratio 
for this case is about the same as that for the case of an uncoated $00 A. 0 
platinum film on pyrex. However, thicker coatings yield larger values for 
the heat flux ratio. 

Variation of substrate properties with temperature: An assumption in the 

analysis leading to Eq. (5S) and then to Eq. (61) was that the substrate 
properties are constant. Since the substrate does undergo a change in temper- 
ature during testing, the accompanying change in substrate property variation 
could have an influence. This assumption has been examined and results are 
reported in Refs. 19 and 21. Fig. 11 summarizes the results for a platinum 
film on pyrex for two boundary conditions in terms of the ratio q^, the heat 
flux computed using variable substrate properties, to q^, the heat flux 
computed assuming constant substrate properties. The boundary conditions 
noted in Fig. 11 are approximately those for q(0,t) = const (boundary condi- 
tion (a)) and heat flux prescribed by q (0, t = const (boundary condition (b)). 
The gage substrate surface temperature change encountered in shock tube flows 
over surfaces parallel to the flow direction such as that over a flat plate 
can be of the order of tens of degrees Centigrade. It can be seen from Fig. 11 
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that the variation of substrate properties should not be neglected. In 
general incorporation of variable substrate properties into a numerical 
procedure as v?a 3 done in the preparation of the curves for Fig. 11 is a 
complex matter (see Refs. 19 and 21). Therefore it is desirable to have at 
hand a simple means for dealing ■with variable properties. A simple but 
accurate means of incorporating variable properties into the determination 
of heat flux involves use of Eq. (6l), the results in Fig. 11, the temperature 
function » an d a fictitious variation of substrate thermal conductivity 

t 

k with temperature. Assuming a linear variation of k vfith T, J3 becomes 

i $(0,t) = T(0,t) Q 1 + (S/2k p ) T(0,t)] (61a) 

where S is the slope of the k vs T curve, S = 0 corresponds to the constant 
substrate property case yielding 0 ~ T(0,t), which is the temperature 
variation with time that is used to compute the variation of heat flux with 
time for the constant properties case using Eq. (6l). However, if S/2k^ is 
taken as 0.0020/°C and if the temperature f(0,t) determined from Eq. (61a) is 
substituted for T(0,t) in Eq. (6l), the resulting heat flux values when divided 
by the corresponding .values computed assuming constant substrate properties 
will yield a curve q v /q^ that is approximately an average of curves (a) and 
(b) in Fig. 11. Hence by use of the fictitious surface temperature change 
£j(0,t) in Eq. (6^ the effects of variable substrate properties on heat flux 
determination can be reasonably accounted for. 
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Summary of factors influencing heat flux determination . — Since several 
factors have a bearing on evaluation of heat flux from gage response it is 
of interest to present a summary with comments on the relative importance 
of each. 

*ff Definition of response curve: First and foremost it is essential that 

the gage response be recorded and displayed in such a manner that the 
response be clearly defined. The horizontal and vertical sensitivi- 
ties of the oscilloscope should be selected to obtain as large a 
trace as possible, consistent with the usable range of the oscillo- 
•'scope. Also the Polaroid photo of the trace and the oscilloscope grid 
should be carefully taken. If early as well as lace time evaluation 
of the heat flux is to be made two oscilloscopes in parallel should be 
used to display the gage response, one with a faster sweep time to 
accurately define the early time portion of the response curve. Two 
methods are available that work well for determining discrete values 
from the response curve for heat flux evaluation. Use of an optical 
reader which enlarges the trace and permits coordinates to be read 
works well when available. A second satisfactory method consists 
of making a transparent copy of the Polaroid photo of the gage 
response and forming an enlarged image of the photo on graph paper 
by means of a projector. The trace and grid are then transferred to 
the graph paper and values that define the response curve are read 
from the graph paper. 



Calibration: Calibration quantities related to the gage are P in 

Eq. (57) and \J ( k p c ) s in Eq. (61). Accurate determination of 
these quantities is important since the heat flux is inversely propor- 
tional and proportional respectively to these quantities. Experience 
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has shown that it is necessary to determine P for each individual gage. 

This is done by measuring resistance vs temperature for the gage using a 
constant temperature bath. As previously mentioned, the accepted value for' 

Y(h c) G is 0.0737 Btu/ft^ °F sec' 2 . Experience has shown that individual 
calibration of each gage to determine \/(k p c) g is unnecessary. 

Thin- film and coating thickness 1 After care has been taken to obtain good 
trace definition and proper calibration, corrections for film and coating 
thicknesses should be considered if very early times are of interest or if 
relative thick coatings or films are used. Although the combined influence 
of thickness of the film and coating together has not been analyzed, inferences 
can be made from the analysis of the individual effects « It is noted that 
the heat flux ratio curve for a platinum film on pyrex in Fig, 9(b) is for a 
film thickness of ^>00 A , which is the approximate upper limit expected for fiXm 
thickness. Therefore, the correction for platinum film thickness for most 
gages would be less than those inferred from Fig, 9(b) and would be negligible 
in practice except at very early times. However, when a coating is used, 
the coating thickness may not be negligible. This can be noted from the fact 
that the 1000 A° coating thickness assumed for Fig. 9(b) (pyrex on pyrex) 
represents the approximate lower limit of coating thickness necessary to 
protect the platinum film from ionization in the gas; and, as previously 
noted, coatings with thicknesses greater than 1000 A 0 yield heat flux ratio 
curves that are higher than the one in Fig, 9(b), This indicates that 
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significant errors may be introduced if corrections for the coating thickness 

(like those made using Fig. 8 for the example of Fig. 9) are not included. 

When corrections for film or coating thickness are applied, it is more 

important' that the thickness be known than it is to know the nature of the 

< 

surface heat flux variation with time. 

Variable substrate properties: This influence is important only for large 

changes in surface temperature as noted in Fig. 11. However, as previously 
noted, the influence of variable substrate properties can readily be incorpor 
ated by use of the temperature function $(0,t) from Eq. (61a) (with S/2k^ = 
0.0018/°C) in place of T(0,t) in Eq. (6l). 

An accurate assessment of the overall uncertainty involved in measuring 

heat transfer rates using thin film gages is difficult to make due to the 

large number of factors (and the uncertainty within these factors) that 

influence the results. However, estimates made on the basis of the method 

25 

Kline and McClintock and experience indicate that with reasonable care heat 


flux can be measured within +8%. 
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EXPERIMENTAL INVESTIGATION 
Shock Tubes 

Three shock tubes were employed in obtaining experimental data in this 

investigation. The important quantities re.lated to the use of the shock 

tube in this study are the shock speed U and the initial temperature and 

s 

pressure of the test gas. In all cases the initial test gas temperature v/as 
assumed to be room temperature. The following is a brief description of 
the three shock tubes used. 

Iowa State Engineering Research Institute Shock Tube . — This shock 
tube is cold-gas driven and has a 3 in. x 6 in. rectangular cross section. 

The test model was located 30 feet from the diaphragm. Shock speeds 
are measured using microsecond interval counters which are triggered 
by the amplified output of thin- film transducers mounted flush v?ith the tube 
side wall and positioned at 1 foot intervals just upstream of the test 
section. For the data taken in this investigation the initial test-gas 
pressure was measured using a Veco thermocouple-type vacuum gage which was 
calibrated using McLeod gages. 

Three- Inch Combustion-driven Shock Tube . — This shock tube is located 
at the NASA-Ames Research Center. It can be operated in either the cold-gas 
or combustion driven mode and has a circular cross section three inches in 
diameter. The test model was located 25 feet from the diaphragm. The shock 
speed was measured using microsecond interval counters which were triggered 
by the amplified output of pressure transducers spaced at known distances 
in the shock tube side wall. The initial test gas pressure p^ was measured 
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by means of a Wallace and Tiernan bourdon-type absolute pressure gage in the 
range of about 5 torr and a McLeod gage for pressures in the range of about 
one torr. 

Six-inch Arc-driven Shock Tube . --This shock tube is also located at 
the NASA-Ames Research Center. The test model was located 25 ft. from the 
diaphragm. Shock speeds were measured by means of ionization probes, associ- 
ated electronics, and a microsecond counter. The probes were spaced 2 feet 
apart at a distance 1.5 upstream of the model. The desired initial test gas 
pressure p^ was obtained using the gas-loading device which is an integral 
part of the shock tube facility. The gas-loading device was calibrated using 
a McLeod gage. 
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Models 

Flat plate models were designed to completely span the shock tube 
except for a small clearance and were mounted in the tube with zero angle 
of attack with the plate test surface approximately on the tube centerline. 
Fig. 12 schematically displays the flat plate model design. The plate leading 
edge angle to in Fig. 12 was chosen small enough to attach the leading edge 
shock. The leading edge radius was maintained at approximately 0.002 inches 
in order to minimize the disturbance of the flow over the plate upper surface 
caused by leading edge bluntness. The plate length was chosen consistent 
with the expected hot-flow testing times and the desire to obtain data over a 
wide range of the distance- time variable a. The important features of the 
plate models are listed in Table I. 

Thin-film heat flux gages (see section on Heat Transfer Measurement) 

were positioned at two or more x locations, x , on the plate. These gages 

g 

were in most cases prepared using substrate elements approximately l/2 ,r x 1” x 
1/8” thick. Platinum films were applied by the painting and firing technique 
using Hanovia 05-x Platinum Bright. Three coats of Platinum Bright were 
applied to each gage. Table II lists the gages that were used to obtain the 
heat transfer rates presented in this report. Care was taken to mount the 
gages flush with the plate test surface in the manner noted in Fig. 12. 

Potting materials were used to fill any major surface imperfections in the 
region bordering the gages. In some tests heat flux gages were also mounted 
flush with the shock tube side wall to measure the heat flux through the 
side-wall boundary layer. This measurement was made to provide additional 
comparative data to verify plate heat transfer measurements in the region of 
the flat plate boundary layer (the M region) that is identical to that on 


the side wall. 
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TABLE 1. FIAT PLATE MODELS 


Plate 

No. 

Width x 
Length 

Gage 

Locations, x.. 

co ; 

Fig 12(b) 

Shock 

Tube 

Plate Mounting 
Method 

1 

3 in. x 
6 in. 

* O 

35.5, 60.5, 

86 . 5 mm 

14° 

NASA 3 in. 

I.S.U. 

cantilevered 
from end wall 

2 

3 in. x 
6 in. 

52,77 mm 

14° 

NASA 3 in. 

cantilevered 
from end wall 

3 

6 in . x 
8 ini. 

Set; k; 14.7 
28.4, 75.4 mm 
Set B ; 15.1 
58.6, 89 . 9 trim 

25° 

NASA 6 in. 

anchored by 
pins through 
side wall 

1 
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TABLE II. THIN FILM HEAT FLUX GAGES 


Gage 
No . 

Resistance R 
ohms at 74°F° 

r,Eq.(57), 
°R" 1 

Coating 

Use 

2 

56.3 

0.00141 

None 

Plate 2, Xg » 52 mm 

3 

74.6 

0.00142 

None 

Plate 2 , x = 77 mm 
8 

15 

111.6 

0.00124 

None 

Plate 1 , Xg *» 35 . 5 mm 

16 

97.9 

0.00131 

None 

Plate 1, x = 60.5 rnm 

o 

17 

116.8 

0.00125 

None 

Plate 1, x *= 86.5 mm 

O 

25 

90.1 

0.00125 

SiO, 600A° 

Sidewall, ISU Shock Tube 

40 

124.8 

0.00122 

None 

Sidewall, ISU Shock Tube 

M4 

106 .2 

0.00135 

a SiO, 2000A° 

Plate 3, x = 14.7 ran 
g 

N3 

203.9 

0.00134 

a SiO, 4000A° 

Plate 3, x = 58.6 mm 

O 

N4 

157.2 

0.00137 

a SiO, 4000A° 

Plate 3, x " 88.9 mm 
& 

T3 

142.0 

0.00130 

c Sapphire 

Plate 3, x = 28.4 mm 

O 

T9 

191.6 

0.00135 

a Si0, 2000A° 

Plate 3, x - 74.5 mnt 
s 

V9 

143.4 

0.00133 

a SiO, 4000A° 

Plate 3, x » 15.1 mm 
g 

02 

106.3 

0.00140 

b SiO, 900A° 

Plate 1, x 35.5 mm 

O 

V? 

160.0 

0.00135 

c 5apphire 

Plate 1, x » 35.5 nun 

V8 

155.0 

0.00135 

c Sapphire 

! 

Plate 1, Xg = 60.5ro:n 


a Placed in oven for 16 hrs. at 900°F to oxidize SiO to Si02. 
^Placed in oven for 6 hrs. at 900°F to oxidize SiO to Si02> . 


. c Thickness unknown 
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Data Analysis 

For a given shock tube test the following dataware obtained i 

1. The initial pressure and temperature, T l» £or the 

test gas. 

2. The shock speed U g . * 

3. Oscilloscope records of the response of each thin-film heat- 
flux gage. 

Using these data and appropriate calibration quantities, experimental 
values of Q vs o< were computed in order to make comparisons with the 
corresponding theoretical values of Q y vs c< (Fig* 3) determined foi the 
particular shock tube test. Experimental values of Q x vr> were computed 

by Eq. (53) which may be written as 



The sources of the various quantities in Eq. (62) are discussed below. 

Heat transfer rates q w =■ q w (t') at each gage location x = Xg were 
determined from gage response recorded by oscilloscopes, using the numerical 
method of Eq. (61) and the response curve definition techniques described 
under Definition of resp o nse curve s in the section of this report dealing 
with heat flux measurement. Note that t in Eq. (61) must be interpreted as 
t r , where t' is zero when the shock arrives at the gage. The resistance- 
temperature quantity F in Eq. (57) v?as determined for each heat flux gage 
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using measured values of gage resistance vs temperature. The variation in 
gage substrate properties was accounted for using the previously-described 
method embodied in Eq, (61a). 

Corrections for the platinum film thickness were not made to gage response 
curves since as previously noted the error involved in ignoring the film 
thickness appears to be quite small. However, as noted in the discussion 
of Fig. 9b, neglecting the presence of the coatings might lead to larger 
errors. Therefore the presence of the coating was not neglected. The 
influence of the coatings on gage response could not be accounted for in an 
exact manner for two reasons. First, only the approximate thickness of the 
coatings (as indicated by the thickness monitor on the vacuum deposition 
apparatus) v?ere known. A profilometer was not available to measure the 
coating thickness. Second, the nature of most of the coatings used was not 
known with certainty. As noted in Table II, two materials were used in 
coating the gages: SiO and sapphire. Mo9t of the results presented in this 
report where coated gages were used were obtained using gages coated with SiO. 
Although sapphire has desirable thermal characteristics and appears to offer 
good protection against ionization effects in the flow, it was not used 
extensively due to difficulties encountered in obtaining bonding of the 
sapphire coatings to the gage. The gages with SiO coatings were 

placed in an oven to oxidize the SiO to Si0 2 . Evidence exists (Ref. 26) 
that SiO coatings in excess of 1000 A° in thickness never completely oxidize 
due to the formation of a crust of Si0 2 that inhibits oxidization of the 
remainder of the coating. Therefore due to. the uncertainties in both thick.-. , 
nesses and the nature of the SiO coatings and the absence of other information 
the following assumptions were made: the thicknesses listed in Table II 
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were correct, and the thermal properties of the coatings were those of 
pyrex. Corrections for the sapphire coatings were not made in the two 
instances where such coatings were used, because, due to the desirable 
thermal characteristics of the sapphire, corrections appear to be unnecessary. 
With these assumptions, values for the coating correction factor cf for Eq. (59) 
were determined from curves with ^ " 1 and qVt" = const (like those in Fig. 8) 
plotted for the given coating thickness. This correction coupled with the 
correction associated with variable properties for pyrex constituted 
the two corrections that were incorporated into the determination of heat 
transfer rates. A further assessment of the influence of these corrections 
on experimental results is made later in this report. 

The values p^ and T^ in Eq. (62) are the measured values for the initial 
pressure and temperature of the test gas and » U^/A, where A “‘ i )* 

The quantity C w » $ w }** x/ f oj^o determined from real-gas 

thermo dynamic and transport property relations for the test gas. The ratio 
h /h^ is the enthalpy ratio across the shock and is determined from the 
appropriate real gas normal shock chart using the measured values for U^, 
and T^. JXy is the viscosity of the test gas evaluated at T^, and is 

the Prandtl number at T w ~ f> Q f was also determined from the 

appropriate normal shock chart. 

Values of corresponding to the values determined from Eq. (&2) are 


y _ ^3 A 

u. 0 t = U 5 U' + V^) 


*3 A 
h 5 t'ty 3 


( 63 ) 


Thus Eqs. (62) and (63) yield and corresponding o: values at each gage 
location that can be compared with theoretical predictions for vs oC . 
As noted previously, the theory for the boundary layer under consideration 
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predicts a single curve Q vs ex' for all x locations on the plate. Therefore 
the gages located at several x^ values provide data for verification of 
this key aspect of the theory. 



47 


RESULTS AND DISCUSSION 

Shock tube experiments were carried out in this investigation that 
produced experimental heat transfer data for g 5 flows induced by shocks 
with speeds ranging from 1.695 to 7.34 km/ sec. .Air and nitrogen were used 

f 

as test gases. Experiments with air were carried out for shock speeds 
ranging from 1.695 to 2.83 km/sec. All heat transfer measurements made in 
air were made with uncoated thin-film heat-flux gages. At shock speeds 
greater than about 3 km/sec in air or nitrogen, ionization of the test gas 

I 

! 

induces spurious gage response signals when uncoated gages are used. Experi- 
ments with nitrogen as the test gas were carried out using coated gages for 
shock speeds ranging from 3 <16 to 7*34 km/ sec. The results obtained for 
each test gas will be presented and discussed separately* 

Experiments in Air (Uncoated Heat-flux Gages) 

The results of the experiments with air as the test gas are compared 

with theoretical predictions in Figs. 13 to 22. These tests were made using 

the Iowa State chock tube and the NASA 3 inch shock tube. The theoretical 

curves for vs ei in these figures were determined using the ratios from 

Fig. 4 and Eqs. (54) and (55). The steady state value, Eq. (32), and 

at oC « 1 for Eq. (34) were determined from Fig. 5 for shock speeds below 2 

km/sec. For higher shock speeds these quantities ware obtained by solving 

Eqs. (22) and (23) for each case using the Praivdtl number variation from Ref. 

8 and the fl* variation from Ref. 9. This was necessary because the pressure 

behind the shock at the higher shock speeds differed considerably from the 

pressure of one atmosphere to which Fig. 5. is limited. 

The theoretical curves Q , vs ■oi in Figs. 13 to 22 are strictly applicable 

only for the condition U t £ L, i.e., the flow- initiating shock wave (Fig. la) 

s 
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remains on the plate. For a given heat flux gage, the value of t' at which 

the shock reaches the trailing edge of the plate is t£ ~ (L - x^) /U g . (See 

Fig. 12) The corresponding value can be found from Eq. (63). The 

condition U t =* L is denoted by the half-filled symbol in Figs. 13 to 22. 
s i 

A heat flux gage response record for* the tests made in air is shown in 
Fig. 23a. This record, which was obtained in Run C 142 for gage 15 (Fig. 22) 
at a shock speed of 2.83 km/sec^ is typical of those obtained for the air 
tests and shows no significant influence from ionization in the flow. 

For runs made at appreciably higher shock speeds, ionization of the 

air rendered the gage response signals useless. The region 2 testing times 

as determined by examining the gage response records were in most cases 

appreciably longer than t' for the gage nearest the leading edge of the 

- - - 

plate. 

Fig. 13 presents results for a shock speed of 1.695 km/sec. The extended 
theoretical steady state and M region solutions are shown for purposes 
of comparison with the theoretical interaction region solution. The experi- 
mental results generally agree well with the theory at lower values of c4 , 
but tend to disagree with the theory at larger values of ©£, Fig. 14(a) shows 
a similar behavior for comparison of experimental results and theory’ for a 

case with U = 1.966 km/sec. Examination of Eq. (63) indicates that for a 
s 

\ 

given case, small values of t' correspond to large values of Thus the 

region of disagreement in Figs. 13 and 14(a) corresponds to small t'. It is 
therefore of interest to examine the heat transfer results in terms of t 1 . 

Fig. 14(b) presents dimensional heat transfer results for Run N010, Fig. 14(a), 

Simultaneous to measurement of the plate heat 


as a function of t' . 
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transfer presented in Fig. 14(a), the side wall heat flux was measured in 

Run N 010. As indicated by Eq. (30), q w (t') 1 / 2 = const for the side wall 

boundary layer. The theoretical value of q (t*) 1 / 2 for Run N 010 is indicated 

in Fig. 14(b). Also shown are experimental values of q w (t') 1//2 determined from 

the response of thin film heat flux gage 40 which was mounted flush with the 

side wall of the shock tube. Agreement of these values with theory is very good 

after the first few microseconds. Experimental values of q w (t') 1//2 as indicated 

by plate gages 15 and 16 for Run N 010 are also shown in Fig, 14(b). These 

values are of interest because, according to theory, the heat flux measured 

by the plate gages should be identical to that for the side wall boundary 

layer for small values of time after the shock passes over the plate gages, i.e., 

small values of t r . At t' = (x /U )/(A-l), i.e., at a = 1, the side wall region 

g s 

(the M region) terminates and the interaction region begins. For somewhat larger 

values of t ' (values of a less than unity) the theoretical solution is 

essentially that for the extended M region boundai'y layer. With still larger 

values of t* (smaller a) the solution for the heat flux should depart from the 

value of q,,(t') ly/2 for the side wall. This behavior is exhibited for t* 

greater than about 15 usee in Fig. 14(b) by the experimental heat transfer 

rates indicated by plate gages 15 and 16. For a short time interval after 

the first few microseconds the two plate gages indicate values of q (t 1 ) 1/2 

w ■ 

that are in very good agreement with the side wall theory. With increasing 
t’, q w (t') l/2 values from gage 15 (the gage nearest the plate leading edge) 
depart as expected from the side wall theory and are followed by those 
indicated by gage 16. 
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Values of t' at which the M region terminates are indicated in Fig. 

14(b) for both plate gages. It is seen from both Figs. 14(a) and 14(b) that 
the M regions terminate prior to the time t' at which experimental and 

1 t 

theoretical M region values for q^(t')’' 5 agree. For small values of t’ 

(less than about lO^sec), the experimental values of q w ( t ' ) ^ indicated by the 
three gages are in good agreement but are larger than the theory predicts. 

The reason for this behavior is related in part to early time errors in the 
numerical method of Eq. (61) and to difficulties' encountered in defining 
both the xero for t* and the shape of the gage response trace at small values 
of t’. It should be noted that errors of this nature tend to rapidly 
diminish as time increases and therefore do not significantly influence the 
experimental results at larger values of time. 

Fig, 14(b) also presents for Run K 010 a comparison of the theoretical 
steady state value for q w (x)'*, which according to Eq. (53) is constant, and 
experimental values for q w (x)" 2 determined from the heat flux indicated by 
gages 15 and 16. The experimental results for both gages approach the steady 

j,. 

state theoretical value of q w (x) 2 as time increases and the experimental 
values agree well with the theory at steady state. This comparison of results 
as steady state is approached is interesting but not as universal as the 
comparison in terms of Q x and in Fig. 14(a). 

Figs. 15(a) and 15(b) present results for a shock speed of 2.35 km/sec 
in a manner similar to those in Fig. 14(a) and 14(b). The results as presented 
in Fig. 15(a) are in good agreement with the theory. As noted in Fig. 15(b), 
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a side walx gage- vas employed to measure the side wall heat flux. Again, 
except for the first few microseconds, good agreement exists between the 
theoretical M region value for q w (t') 2 and the experimental values for that 
quantity indicated by each of the three gages. For the time 'interval shown in 
Fig. lj(b) only the heat flux quantity for plate gage 15, the gage nearest 
the leading edge, indicates an influence of the interaction region. 

Figures 16 to 22 present results for the remainder of the shock speed 
range covered in this investigation using air as the test gas. In these figures 
experimental results are presented only for values of corresponding to t* 

values larger than the times at which experimental values of C"t f 2 con- 

stant with t f . With the exception of Fig. 16, 18(a) and 20, the experimental 
results are observed to be in good agreement with the theory. In Figs. 16, 18(a) 
and 20 the experimental results generally fall somewhat above the theoretical 
curves. It is believed that this disagreement is due to errors associated with 
the experimental phase of the investigation rather than a failure of the theory 
to properly predict variation of Q x with <X . Possible sources of experimental 
error are numerous, but those believed to be involved here are errors in voltage 
measurement and those involved with measuring the initial pressure p_j_ of the 
test gas. The initial pressure of the test gas appears as 1/^)^ in Eq. (62) 
which was used to compute experimental values of Q . Vertical scale calibration 

of the oscilloscope affects the voltage quantities in Eq. (57) that relate to 

% 

the determination of the gage substrate surface temperature change T(0,t) } 
which in turn is used in the computation of heat flux by Eq. (61). Errors 



52 


ill voltage measurement are directly related to errors in heat flux, i.e., 
a 10% error in. voltage measurement results in a corresponding 10% error in 
heat flux. Evidence that a fixed experimental error is involved with the 
disagreement in Figs. 16, 18(a), and 20 is obtained from the results for 
Run C 137, Fig. 18(a), as presented in Fig. 18(b). The results in this 
figure were obtained by multiplying the corresponding experimental 0^ 
quantities in Fig. 18(a) by the ratio of the theoretical to experimental 
Q x , the ratio being determined at low values of a for each gage. This 
results in good agreement of experimental results with theory at larger 
values of a as can be seen in Fig. 18(b). The results in Figs. 16 and 20 
exhibit the same trends as those in Fig. 18(a) and when modified in a manner 
similar to those in Fig. 18(b) also show good agreement with theory thereby 
providing additional evidence that fixed experimental errors are involved in 
the results of Figs. 16, 18(a), and 20. 

The comparison of results for shock speeds ranging from 2.35 to 2.83 
km/sec in Figs. 15 to 22 are for higher shock speeds than those obtained by 
Felderman^*^ using air as the test gas. 
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Experiments in Nitrogen (Coated Heat-flux Gages) 

The results of experiments performed with nitrogen as the test gas are 

compared with theoretical heat transfer predictions in Figs. 24 to 32. These 

r 1 

experiments cover a range of shock speed from 3.16 to 7.34 km/sec, and were 
carried out using two shock tubes, the 3 inch NASA shock tube (combustion 
driven) and the NASA 6 inch arc-driven shock tube. Figures 24 and 25 present 
results obtained using the 3 inch shock tube while Figs. 26 to 32 report the 
results obtained for the 6 inch shock tube. In two instances more than one 
run was made at essentially the same shock speed. These cases are reported 
in Figs. 30(a) and (b) and in Figs. 31(a), (b) and (c) . Results in the 
form q(x)' 5 vs t’ are shown for two selected cases in Figs. 28(b) and 32(b). 

The theoretical curves vs <s4 in Figs, 24 to 32 were obtained in a 
manner similar to that previously discussed for air in Figs. 13 to 22. The 
steady state and N region values for Q, (Eqs. (32) and (34) respectively) 

A 

were first obtained by numerically solving Eqs. (22) and (23) for both the 
steady state and M regions. The property ratios across the incident shock 
and the variation of properties with enthalpy for nitrogen required for these 
solutions were determined at the appropriate pressure from the charts and 
tables for equilibrium nitrogen presented by Ahtye and Peng, Ref. 27. A 
given curve Q vs<s( was then constructed using the method related to Fig. 3 
and Eqs. (54) and (55). The ratios R gg and R^ were obtained from Fig. 4. 

Fig. 4 is strictly applicable only to air at p “ 1 atm. However, the results 
obtained using this figure for nitrogen should be acceptable first, because 
the denominators of Eqs. (54) and (55) were accurately determined for nitrogen, 
and second, due to the similarities between nitrogen and air, the ratios R 
and R for nitrogen would not be expected to differ significantly from those 
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for air. Use of Fig. 4 for nitrogen eliminated the lengthy and expensive 
procedure required to solve the governing equations for the unsteady 
interaction region, Eqs. (49) and (50), which lead to Eq. (52). 

Figures 23b to 23f show representative heat flux gage response records 
for the tests in nitrogen. Gage coatings for the gages noted are described 
in Table II. All records in this group except Fig. 22b show some influence 
of ionization of the nitrogen on the gage signal. This influence is evident 
in two respects. First, precursory ionization influences the signal prior 
to arrival of the shock at the gage as evidenced by the down-going portion 

I 

of the trace prior to shock arrival in Figs. 23 d, e, and f to a smaller 

i 

extent in Fig. 23c. Second, there is some influence of ionization immediately 

i 

after shock arrival, and to a lesser extent at later times. These effects 
are particularly evident in Fig. 23f in which the effects of precursory 
ionization prior to shock arrival and the somewhat eratic signal after shock 
arrival both exist to a significant extent. 

Figures 23b to 23f illustrate the degree to which various oxidized SiO 
coatings provide protection against ionization. The coating for gage 02 
(900 A° SiO oxidized 6 hrs. at 900°F) provided excellent protection at the 
flow conditions noted in -Fig. 23b. However, at appreciably higher shock speeds eratic 
signals resulted from gage 02. Gage V9, 4000 A 0 SiO coatings oxidized 16 
hrs. at 900°F, gave relatively noise-free response (Figs. 23 c, d, and e) as 
compared to the response of gage M4, Fig. 23f, which was prepared in a 
similar manner except the thickness of the SiO was 2000 A°. The flows related 
to Figs, 23e and 23f involved about the same level of ionization since the 
shock speed and pressure level was about the same. The response curves of 
gage V9 are typical of the response curves of the gages that were used to 
obtain most of the heat flux data for nitrogen flows. Due to the ionization 
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effects on gage response in the immediate vicinity of the shock, heat 
transfer results at very early {rimes after shock arrival at the gage were 
not considered to be reliable. 

Termination of region 2 shock tube flow is in many cases evident from 
the response of heat flux gages. Figures 23d, e, and f clearly show this. 

The eratic signal recorded just prior to the smooth declining segment of 
the response curve is interpreted to be due to the arrival cf the contact 
front that forms between tho test and driver gases. In view of the relatively 
short testing times associated with high shock speeds, a comparison of the 
region 2 testing time with the time required for the shock to reach the trail- 
ing edge of the plate is of interest. Figure 33 presents such a comparison 
for the range of shock speeds covered in tests with Plate 3 (see Table II) 
mounted in the 6 inch arc-driven shock tube. The testing times were taken 
as the time between shock arrival and interface arrival at the front gage on 
the plate. Values of t^ => (L-Xg)/U g for the front gage (Xg ** 15.1 mm) are 
also shown in the Figure. Since the testing time is somewhat greater than 
for the full range of shock speed shown it was concluded that region 2 

was of greater length than the plate and that times up to at least t' for 

L 

each gage involved region 2 flow over the whole plate. 

With a few exceptions the experimental heat transfer results presented in 
Figs. 24 to 32 exhibit good agreement with theoretical predictions for the 
range of shown. Results are not presented in some cases for values of 
corresponding to t' values less than about 5 yjsec because of uncertainty 
regarding the gage response configuration at early times due to ionization 
effects in the vicinity of shock wave arrival. Small values of t' correspond 
to large values of ,c>d (see Eq. <63)). Accordingly, this eliminated some 
data at large values of for each gage location. Of the twelve comparisons 
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of vs o 4 , Figs. 26, 30(a), and 32(a) tend to exhibit experimental 
results somewhat above the theory, while in Figs, 31(a), (b) , and (c) the 
results tend to fall below the theory. In the remaining six figures the 
experimental results generally agree well with the theory. Agreement between 
the results indicated by different gages varies from quite good in Fig. 28(a) 
to fair in Fig. 31(b). In Figs. 26 to 32 Gage V9 (the front gage) tends to 
yield results a few percent below the other two gages. Figure 34 presents a 

comparison of experimental steady state values of with theoretical steady 

! ' • 

state values for the range of shock speeds covered using nitrogen as the 

test gas. In this figure only experimental results for those gages 
indicating a steady heat flux with time are shown. These results agree with 
theory within about + 107». 

As previously noted the theoretical curves vs C>4 are strictly 
applicable only as long as the shock remains on the plate, i.e., t ^ k/U s , 
and correspondingly, t' ^ (L - x^)/U g . The half-filled symbols for the experi- 
mental results for nitrogen in Figs. 24 to 32 and for air in Figs. 13 to 22 
indicate the condition t « L/U s . In many of these figures this condition 
occurs at values of c*£ in the unsteady interaction region. Experimental 
results at values of o<£ smaller than those corresponding to the half-filled sym- 
bols are then related to conditions that exist after the shock has reached the 
end of the plate. Experimental results at these lower values of ©< are shown 
in the figures down to values of o< at which a significant departure from 
the trend of the results occurred. In many cases this departure corresponded 
to the termination of the available test time. Thus the theory for the 
unsteady laminar boundary layer presented here continues to predict the plate 
heat transfer for times significantly longer than t a L/U s . 
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Figures 28b and 32b present results in the form q t vs t' for two 
typical shock speeds; one at an intermediate speed of 4.59 km/sec (fig* 28b) 
and the other at the highest shock speed investigated, 7,34 km/sec (Fig- 32b). 
Experimental results are shown only for those heat flux gages which indicated 
steady state heat flux values. The apparent "noise" is due mostly to the 
date reduction technique. Typical influence of the two corrections applied 
in reducing the thin film gage response data is illustrated in these figures 
for gage V9. The correction for the coating has the largest influence at 
small values of t', while the variable substrate properties correction is 
related to temperature through t' and increases with t' since surface tempera- 
ture increases with t 1 . Both corrections are of the order of a few percent 
in Fig. 28b. Ignoring the corrections would lower the results by 8 to 10 
percent. In Fig. 32b the variable substrate properties correction is larger 
due to the larger change in substrate surface temperature accompanying the 
higher heat flux encountered for the higher shock speed. Ignoring the 
corrections here would yield results about 15% lower. Although both of the 
corrections made are approximate, they seem to be justified by the experimental 


results . 
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CONCLUDING REMARKS 

The comparisons of theoretical and experimental heat transfer rates 

made for the range of shock speeds covered in this study generally verify 

the theory of the shock- induced unsteady laminar boundary layer on a flat 

* 

plate as presented here and in Ref. 1. In particular, the results confirm 
the self similarity that is revealed by the theoretical analysis, namely 
that the wall heat transfer and other boundary layer quantities can be 
formulated in terms of the time-position variable alone. The theory might 
be expected to predict the boundary layer behavior only for the time interval 
during which incident shock remains on the plate. However, the experimental 
results indicate that the theory continues to predict the boundary layer 
behavior for significantly longer intervals of time. 

Although some difficulties were encountered in measuring heat transfer 
rates at values of oC in the vicinity of unity, the heat transfer measurement 
technique employed was generally satisfactory* In the cases where measured 
heat transfer rates differed to a pronounced degree with the theory, there 
is evidence that the difference was due to some fixed error in an experi- 
mentally measured quantity rather than in a gross failure of the theory to 
predict the heat transfer rates. Some refinements are needed in the heat 
transfer measurement technique, particularly in those related to measurement 

of heat flux at higher shock speeds where test gas ionization influences the 

♦ 

performance of thin- film heat flux gages. Development of gage coatings to 
provide protection against high levels of ionization and interpretation of 
coated gage response to determine heat transfer rates should be the subject 
of further investigation. 

Finally, the experimental results of this investigation confirm that 
under certain conditions steady flow can readily be attained over models 



59 


mounted in shock tubes, even for short-duration flows generated by high 

speed shocks. This fact coupled with' the technique described here for 

measuring time dependent heat flux makes feasible experimental heat transfer 

studies for both laminar and turbulent boundary layer flows over models 

< 

under conditions of high wall cooling and low supersonic free-stream Mach 


numbers . 
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(a) The shock-induced laminar boundary layer on a flat plate. 



(b) The coordinates of Lam and Crocco for the unsteady shock-induced 
laminar boundary layer on a flat plate. 

Figure 1.- Nomenclature for the shock-induced laminar boundary layer on 

a flat plate. 







(b) M region boundary layer in wall-fixed coordinates. 
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Figure 3*- Typical solution for vs a. 
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Figure 7,- Dimensionless interface temperature change of a composite semi- 

infinite solid* 





Assumed gage temperature response 
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from eq. (6l) for 5p = 5(j = 0 


Time increment m 
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(a) Assumed gage temperature response and corresponding heat flux quantities 


•Uncoated 500 A° platinum film on pyrex substrate 

y—1000 A 0 pyrex coating on zero thickness 
. / platinum film with pyrex substrate 


i-i° -\\ , 

,t)c V/ 

^)gr „ _ n: 


t, |isec 

(b) Ratio of corrected to uncorrected heat flux for assumed gage response 
Figure 9.- Assumed gage response curve and computed heat flux quantities. 




Percent difference from exact heat flux 



Figure 10.- Percent error in heat flux involved in using eq. (6l) vs number 
of equal time intervals for two boundary conditions, (ref. 19), 


Boundary conditions 
(a) o(0,t) = cons t(t l/E ) 
(To ) $(0,t) = const 


Surface temperature change, C 

Figure 11.- Ratio of heat flux computed for variable substrate prop 
erties to that computed assuming constant substrate properties . 
Code 77*1-0 pyrex, (ref. 2l). 
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(b) Side view 

Figure 12.- Flat plate model showing typical thin-film heat flux gage 

installation. 



Figure 13.- Comparison of theoretical and experimental heat transfer results 

Uo = I.695 ton/ sec. 



Run N 010, air 

ISU shock tube, Re/x = 613/mm 
U s = I.966 km/sec, M s = 5*66 

p = 2.38 torr, Ti = 300 °K 
p Q = 0.12 atm, A = 1.202 

P Q /P 1 - 5.95 

Po/Pj. = 38.4 .ref. 11 

ho/hi - T»25 

Plate 1 

Gage 15, Xg = 35.5 mm 
Gage 16, Xg = 60.5 mm 
Theory 

Half filled symbols indicate shock 
arrival at plate trailing edge 




O Gage 15 , plate 
□ Gage 16, plate 
o Gage hO, side wall 


Half filled symbols indicate shock 
arrival at plate trailing edge 


q^x 1 ^ 2 = 9*93 Btu/ft s ^ 2 sec, steady state theory 




■Termination of M region, gage 15 


q (t*) 1 ^ 2 = O.I985 Btu/ft^ec 1 ^ 2 , side wall theory 


Termination of M region, gage 16 

I . I . I 


(b) Heat transfer quantities vs t 


Figure 14.- Concluded. 




Run N 012, air 


ISU shock tube, Re/x = 202/ni 
u s = 2.35 km/sec M s = 6.78 

p x - O .63 torr, Ti = 3 OO °K 
p D = 0.0476 atm, A = 1.172 


P Q p 1 = 6.85 1 

P Q /Pi = 57.3 
faoAi = 10 . 01 . 


ref. 11 


Plate 1 

Gage 15 , Xg = 35*5 nnn 
Gage 16, Xg = 6 O .5 mill 
Theory- 

Half filled symbols indicate 
shock arrival at plate 
trailing edge 


(a) Comparison of theoretical and experimental heat transfer results 
Figure 15 •- Data for U s = 2.35 km/sec. 






Figure 15 . - Concluded. 



0 *2 • ^ < 6 »8 1*0 

a = x/aij-jt 


Figure IT.- Comparison of theoretical and experimental heat transfer results 



Run C 137, air 

RASA 3 inch shock tube, Re/x = 660/ir 
U s = 2.56 km/sec, M s = 7.^ 
p = 2.01 torr, T x = 295 °K 
p D - 0.18 atm, A = 1.158 

P 0 /P 1 = 7.32 

p 0 /p i = 68.1 • ref. 11 

ho/h x = 11.83 . 

Plate 1 

Gage 15, x g = 35.5 ram 
Gage 16, Xg = 60.5 mm 
Gage 17, Xg - 86.5 
Theory 


Half filled symbols indicate 
shock arrival at trailing 
edge 


o 

□ O 
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(a) Comparison of theoretical and experimental heat transfer results 
Figure 18.- Data for U s = 2.56 km/sec. 



Flagged symbols are for run C 134 


Half filled symbols indicate shock 
arrival at plate trailing edge 




(b) Modified experimental results 
Figure 18.- Concluded. 



Run N Ol4, air 

ISU shock tube, Re/x = 132/mm 
U s = 2.61 km/sec, M s = 7*5^ 

Pl = 0.52 torr, T x = 298 °K 
p Q = 0.0482 atm, A = 1.148 

P 0 / Pl = 7.67 

Po / Pi = 70.3 -ref. 11 

h Q /h 1 = 12.16 

Plate 1 

Gage 15, x g - 35.5 mm 
Gage 16, xg = 6O.5 mm 

Theory , 


Half filled symbols indicate shock 
arrival at plate trailing edge 


Figure 19.- Comparison of theoretical and experimental heat transfer 

results, U s = 2.61 km/sec. 



Run C 135, air 


NASA shock tube, Re/x = ’J^b/xm 
U s = 2.63 km/sec, M s = 7*62 
p x - 2.10 torr, = 297 °K 
P 0 “ 0*199 atm, A = I.I35 

P 0 /P x = 7*^5 

vjv^ = 72.0 • ref. 11 

h 0 /h 1 = 12 . bO , 

Plate 1 

Gage 15 , Xg = 35*5 nun 
Gage 16 , Xg = 6O.5 mm 
Gage 17, Xg = 86.5 mm 

Theory- 

Half filled symbols indicate shock 
arrival at plate trailing edge 
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Figure 20 .- Comparison of theoretical and experimental heat transfer 

results. U s = 2.61 km /sec. 





Run C 125, air 

NASA 3 inch shock tube, Re/x = 73®/ 
U s = 2.68 km/sec, M g - 7*76 
= 2.0 torr, T x = 296 °K 
p Q - 0.195 atm, A = 1.150 


P 0 / Pl = 7.65 
P 0 /p x =7^-1 

h Q /h 1 = 12.76 


ref. 11 


Plate 2 

Gage 2, Xg = ^2 mm 
Gage 3, Xg = 77 mm 
Theory 

Half filled symbols indicate shock 
arrival at plate trailing edge 


a = x/v^t 


Figure 21.- Comparison of theoretical and experimental heat transfer 

results. U,. = 2.68 km/sec. 



Hun C 142, air 


NASA 3 inch shock tube, Re/x ~ 870 / 
U s = 2.83 km/sec, M s = 8.20 
P 1 = 2.25 torr, T x = 295 °K 
p Q = 0.246 atm, A = 1.143 

P 0 / Pl = 8.05 ' 

P 0 /p 1 = 83.1 ' ref. 11 
h 0 /h 1 = 14.2 . 

Plate 1 

Gage 15, scope 1, Xg - 35*5 nun 
Gage 15, scope 2, 

Gage 16, x = 6 O .5 mm 

O 

Theory 

Half filled symbols indicate 
shock arrival at plate trailing 


Figure 22.- Comparison of theoretical and experimental heat transfer 

results. TJ„ = 2.8Tldn/sec . 





(a) Hun C 142 {fig. 22) 

Us - 2.83 km/sec, gage 15 
Test gas: air 
p c =« 0.246 atm 
Time base: 19.2 psec/div. 
Vertical: 5*52°C/div. 


(b) Run C 159 (fig. 24) 

Us = 3* 18 km/sec, gage 02 
Test gas: nitrogen, 
p 0 = 0.203 atra 
Time base: 10 psec/div. 
Vertical: 6.20°c/div. 


(c) Run S 182 (fig.' 28(a), (b)) 
U s = 4.59 kro/sec, gage V9 
Test gas: nitrogen, 
p 0 = 0.29 atm 
Time base: 20.4 psec/div.. 
Vertical: 9.8l°c/div. 


Figure 23.- Typical oscilloscope records of thin-film heat flux gage response; all sensitivities 

are per major scale division. 



(d) Run S 184 (fig. 30(a)} 

U s = 5*81 km/sec, gage V9 
Test gas: nitrogen, 
p 0 = 0.473 atra 
Time base: 20.4 psec/div. 
Vertical: 24.6°C/div. 


(e) Run S 188 (fig. 31(a)) 

U s = 6.04 km/sec, gage V9 
Test gas: nitrogen, • 
p 0 = 0.26 atm 
Time base: 20.4 psec/div. 
Vertical: 24.4°C/div. 


(f) Run S 175 (fig. 31(c)) 

Ug - 6.10 km/sec, gage M4 
Test gas: nitrogen, 

Po = 0-537 atm 

Time base: 20.4 psec/div. 

Vertical: 26.6 6 C/div. 


Figure 23.- Concluded. 




Run C 159 > nitrogen 

NASA 3 inch shock tube, Re/x = 50l/mm 
U s = 3.16 kro/sec, M s « 9.1A 
p = 1.5 torr, l’i = 296 °K 
p Q = 0.203 atm, A - I.165 

0 o /p x = 7. 08 

Po / Pl =102.5 ■ ref. 27 
ho/hi = rM7 . 

Plate 1 


Gage 02, x g = 35 O nan 
Theory 

Half filled symbol indicates shock 
arrival at plate trailing edge 




OSS' 


Figure 24*- Comparison of theoretical and experimental heat transfer 

results y U s = 3*16 km/sec. 




Run C 184, nitrogen 

NASA 3 inch shock tube, Re/x = 55 6/m 
U s = 3*43 km/sec, M s = 10.0 
p 1 = 2.0 torr, T x = 296 °K 
p 0 = 0.242 atm, A = 1.155 

Po/Pi =7.44 

p 0 /p x =121.5 ref. 27 

ho/hi =20.64 

Plate 1 

Gage V7, x g = 35 • 5 mm 
Gage V8, Xg = 60.5 mni 
— Theory 

Half filled symbols indicate shock 
arrival at plate trailing edge 



Run S lol, nitrogen 

NASA 6 inch shock tube, Re/x ** hyk/m 
U s = ^*04 km/sec, M s = 11.58 
p 1 = 1.0 torr, T 1 =* 296 °K 
p Q = 0.222 atm, A = 1.121 

P 0 /»i “ 9- 25 I 

Pq/p^ = 168.0 ref. 27 
ho/hi =27-50 
Plate 3 

Gage V9, Xg = 15.1 mm 
Gage N3, x g = 58.6 mm 
Gage x g = 88.9 mm 

— Theory- 

Half filled symbols Indicate shock 
arrival at plate trailing edge 




Run S 189, nitrogen 

HASA 6 inch shock tube, Re/x = 1079/mm 
U s = A. 325 km/sec, M s = 12.39 
^ = 2.0 torr, T ± = 296 °K 
p Q = 0.51 atm, A = 1.112 

Po/ Pl =9.87 
P o/Pi = 194 ref. 27 
h o/ h i = 31.4 

Plate 3 


0 

Gage V9, 

Xg = 15.1 nun 

D 

Gage N3, 

Xg = 58.6 mm 

O 

Gage N4, 

Xg - 88.9 mm 


Theory 

Half filled symbols indicate shock 
arrival at plate trailing edge 


lx 1*5 


O 
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a. = x/uot 


Figure 27.- Comparison of theoretical and experimental heat transfer results 

U s = 4.325 km/sec. 


Run S lo2 ? nitrogen 

NASA 6 inch shock tube, Re/x =* 625/m 
U s = 4.59 km./ sec, M s = 13.15 
Pj_ = 1.0 torr, Tx = 296 'K 
p 0 = 0.29 atm, A = 1.106 

p 0 /p 1 = 10.50 j: 

Po/ Pl = 220 ref. 27 
h o/ h a =35.25 
Plate 3 

Gage V9, Xg = 15.1 mm 
Gage N3, Xg = 58-6 mm 
> Gage N 1 *-, x g - 88.9 mm 
— Theory 

Half filled symbols indicate shock 
arrival at plate trailing edge 
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(a) Comparison of theoretical and experimental heat transfer results 
Figure 28.- Data for U s ~ 4*59 km/sec. 



O Gage V9 

□ Gage H3 

Gage V9, Variable substrate properties 

correction only 

A Gage V9 j no corrections 




qyx 1 / 2 = 109 Btu/ft 3 / 2 
steady state theory 


vwr* aa %• 

^0/ \ A 

I ^ ^ / End of testing time 
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j I 


40 60 

t’, p.sec 
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(b) Comparison of q^x 1 / 2 values. 
Figure 28.- Concluded. 






Figure 29 Comparison of theoretical and experimental heat transfer results 




0 -2 '■ .6 .8 i.o 


a = x/uot 

(a) U s = 5-8l km/sec 


Figure 30.- Comparison of theoretical and experimental heat transfer result 



Figure 30.- Concluded 





0 .2 .b .6 ,8 1.0 

a - x/u 0 t 


( a) U s =* 6.04 km/sec 


Figure 31.- Comparison of theoretical and experimental heat transfer results. 




Run S 186, nitrogen 

NASA 6 inch shock tube, Re/x = 520 /ram 
U s = 6.06 km/sec, M s = 17.36 
Pl = 0.5 torn, T 1 = 296^ 
p Q = O.258 atm, A = 1.077 

P 0 / Pl *» 13.87 

p 0 / Pl = 392 ref. 27 

h 0 /hi =60.9 I 


Plate 3 

Gage V9, Xg 
Gage N3, xg 
Gage R4, Xg 
Theory 


15.1 mm 
58.6 mm 
88.9 nim 


Half filled symbols indicate shock 
arrival at plate trailing edge 
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Figure 31.- Continued 



Run S 175 , nitrogen 

NASA 6 inch shock tube, Re/x = 1005/mm 
U s - 6.10 km/sec, M s = 17*7 
= 1.0 torr, T x - 296 °K 
p 0 = 0.537 atm, A = 1.075 

Po/Pi “ l4 -3 

p 0 /p x = 409 ref. 27 

h o/ h i = 63.5 

Plate 3 

O Gage M4, Xg = 14.7 ram 

□ Gage T3, x g * 28.4 mm 

o Gage T9, Xg = 74.5 mm 

Theory 

Half filled symbols indicate shock 
arrival at plate trailing edge 


Figure 31 •- Concluded 



Run S l87j nitrogen 

NASA 6 inch shock tube, Re/x = 663 /mm 
U s = 7.34 km/sec, M s = 21.03 
p x = 0.5 torr, Tx = 296 °K 
p 0 = 0.382 atm, A = I.O 67 

Pq/Px = 18.0 
P 0 /Px = 581 ref ♦ 2 ? 

h o/ h ! “ 8 9.0 J 

Plate 3 

t Gage V9, Xg = 15.1 nun 
Gage N3, Xg = 58.6 mm 
► Gage N4, Xg - 88.9 mm 
— Theory 

Half filled symbols indicate shock 
arrival at plate trailing edge 
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(a) Comparison of theoretical and experimental heat transfer result 


Figure 32.- Data for Uo » 7.34 km/sec. 







Time, 



Figure 33- - Comparison of t^ values and region 2 testing time for plate 3 
mounted in the NASA 6 inch arc driven shock tube. 




Exper imental results 

o Front gage 
□ Second gage 


o o 




